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The standard model of cosmology is based on the existence of homogeneous surfaces as the 
background arena for structure formation. Homogeneity underpins both general relativistic and 
modified gravity models and is central to the way in which we interpret observations of the CMB 
and the galaxy distribution. However, homogeneity cannot be directly observed in the galaxy 
distribution or CMB, even with perfect observations, since we observe on the past lightcone and not 
on spatial surfaces. We can directly observe and test for isotropy, but to link this to homogeneity, we 
need to assume the Copernican Principle. First, we discuss the link between isotropic observations 
on the past lightcone and isotropic spacetime geometry: what observations do we need to be isotropic 
in order to deduce spacetime isotropy? Second, we discuss what we can say with the Copernican 
assumption. The most powerful result is based on the CMB: the vanishing of the dipole, quadrupole 
and octupole of the CMB is sufficient to impose homogeneity. Real observations lead to near- 
isotropy on large scales - does this lead to near-homogeneity? There are important partial results, 
and we discuss why this remains a difficult open question. Thus we are currently unable to prove 
homogeneity of the Universe on large-scales, even with the Copernican Principle. However we can 
use observations of the CMB, galaxies and clusters to test homogeneity itself. 

I. INTRODUCTION 

The standard model of the Universe - the LCDM "concordance" model - is homogeneous, with structure formation 
described via perturbations. Given the assumption of homogeneity, and if GR correctly describes gravity, then the 
acceleration of the Universe is driven by dark energy. The homogeneous LCDM model is highly successful - a simple, 
predictive model that is compatible with all observations up to now. However, there is still no satisfactory description 
of the dark energy that is central to this model. This motivates the need to probe the foundations of the model. We 
can probe the assumption that GR holds on cosmological scales, by investigating modified gravity theories and by 
devising consistency tests of GR. This probe is only effective if we assume homogeneity. Alternatively, we can assume 
that GR holds and probe the assumption of homogeneity (see also [H[2])- 

A common misconception is that "homogeneity is obvious from the CMB and the galaxy distribution". In fact, we 
cannot directly observe or test homogeneity - since we observe down the past lightcone, and not on spatial surfaces 
that intersect that lightcone (see Fig. [T]). We only see the CMB on a 2-sphere at one redshift, and galaxy surveys 
give us the galaxy distribution on 2-spheres of constant redshift. There are interesting and important analyses of 
the observed galaxy distribution to probe statistical homogeneity (e.g. OH]), but these effectively assume an FLRW 
background geometry. 

What we can directly test for is isotropy of observations. This then raises an important, but often overlooked 
question: what observational quantities need to be isotropic in order to enforce isotropy of the geometry? This 
question is addressed in Sec. [TTJ To answer this question we need a fully nonlinear analysis, since we cannot assume 
a priori any symmetric background spacetime. For observations of the galaxy distribution, the answer is - we need 
isotropic angular distances, number counts, bulk velocities and lensing. Isotropy of the CMB by contrast does not in 
itself enforce spacetime isotropy. 

In order to link isotropy to homogeneity, we have to assume the Copernican Principle (CP), i.e. that we are not at 
a special position in the Universe. The CP is not observationally based; it is an expression of the intrinsic limitation 
of observations from one spacetime location. We consider in Sec. |III| what can be done with the CP. If we have 
isotropy along one worldline, based on the observed galaxy distribution, then the CP leads to homogeneity. A more 
powerful result is that homogeneity follows if all observers see isotropic angular distances up to third order in redshift. 
The strongest basis for homogeneity comes from the CMB. This is often considered to be obvious - but it is far from 
straightforward to show that homogeneity of the metric follows if all observers see isotropic CMB. The proof requires 
the general nonlinear Einstcin-Liouville equations. Remarkably, it is not necessary to assume full CMB isotropy for 
each observer - it is enough that each observer sees isotropy in the CMB only up to the octupole. This is the most 
powerful observational basis for homogeneity currently known. 

Of course the real CMB is not exactly isotropic, but nearly isotropic. Does it follow from near-isotropy that the 
Universe is nearly homogeneous, i.e. perturbed FLRW? This has so far only been shown with further assumptions on 
the gradients and time derivatives of CMB multipoles. 



FIG. 1: We observe down the past lightcone and therefore cannot directly confirm homogeneity. 

It is important to stress from the outset that there are two fundamental limitations: 

• Isotropy and homogeneity of observables can only be meaningfully defined on large enough scales - and the 
nature of the transition scale is only poorly understood. 

• Isotropy and homogeneity of observables can only be meaningfully defined in an average sense - and the problem 
of how to average in general relativity (and other metric theories that are intrinsically nonlinear) remains 
unsolved. 

These unresolved issues are of crucial importance in cosmology, but they are not discussed here; instead, observations 
are treated as idealized. 



II. WHAT IS THE BASIS FOR ISOTROPY? 

Here we consider the situation when the CP is not adopted. What is the observational basis for showing that 
spacetime is isotropic about the worldline of a single observer? Isotropy is directly observable and the best example 
is the CMB, which is isotropic about us to ~ (after the dipole is interpreted as due to our motion relative to 

the cosmic frame, and removed by a boost). Observations of the galaxy distribution do not have the same precision, 
but there is no evidence for anisotropy. First we look at observations of matter and then of the CMB. 

A. Isotropic matter distribution on the past lightcone of one observer 

The dominant cosmological components - cold dark matter and dark energy - have not been independently ob- 
served. Unlike baryonic matter, the dark components are up to now only manifest via their gravitational effect. The 
distribution of dark matter is mapped by weak lensing surveys. But to relate the measured projected potential on 
the sky at each redshift to the dark matter, we require a specific model, such as a perturbed FLRW model. The dark 
matter 4- velocity is usually assumed to be aligned with that of baryonic matter - but this is also based on a perturbed 
FLRW model. 

This unavoidably means that we must impose a model for these dark components - not merely their physical 
properties, but how they relate spatially to observed matter - in order to determine their distribution via cosmological 
observations. A starting point [T] is to assume that the CDM 4- velocity is the same as the baryonic 4- velocity, and 
that we know the primordial ratio of CDM density to baryonic density, as well as the bias factor that relates the 
concentrations of CDM and baryons in clustered matter 



p c known from pb and = u° := u a 



(1) 
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If there is a modified gravity theory that avoids the need for CDM, then we do not need ([I]) - but we are likely to 
need other assumptions on the extra degrees of freedom that mimic CDM. 

If dark energy is in the form of A, then we need to assume that its value is known from non-cosmological physics: 

A known independently of cosmological observations. (2) 

For quintessence and other more complicated forms of dark energy, we would need to assume how the dark energy 
field is distributed in spacetime - since we are not assuming a FLRW geometry a priori. If there is a modified gravity 
theory that avoids the need for dark energy, then ^ is not needed, but assumptions will likely be necessary on the 
extra degrees of freedom that mimic dark energy. 

Given the assumptions on the dark components, what can we say about spacetime geometry if the matter distri- 
bution is isotropic on the past lightcone of the observer? Which observables need to be isotropic in order to deduce 
isotropic geometry? It turns out that 4 independent observables on the lightcone are exactly enough to impose 
isotropy of spacetime. The original result for a baryonic universe [HrE] may be updated to include CDM and A pQ, 
incorporated via |T]) and 

Isotropy of matter distribution on the lightcone — » isotropy of spacetime geometry 

If one fundamental observer comoving with the matter measures isotropy of (a) angular diameter distances, 
(b) number counts, (c) bulk velocities, and (d) lensing, in an expanding dust universe with A, then the 
spacetime is isotropic about the observer's worldline. 

Note that isotropy of bulk velocities is equivalent to vanishing transverse velocities (proper motions) on the observer's 
sky. Isotropy of lensing means that there is no distortion of images, only magnification. 

The proof of this result requires a non-perturbative approach - there is no background to perturb around. Since 
the data is given on the past lightcone of the observer, we need the full general metric, adapted to the past lightcones 
of the observer worldline C. We define observational coordinates x M = (w, y, 9, </>), where x p = (9, </>) are the celestial 
coordinates, w = const are the past light cones on C (y = 0), normalized so that w measures proper time along C, and 
y measures distance down the light rays (w,6,(j>) = const (see Fig. [2]). A convenient choice for y is y — z (redshift) 
on the lightcone of here-and-now, w = Wo, and then keep y comoving with matter off the initial lightcone, so that 
u v = 0. Then the matter 4-velocity and the photon wave-vector are 

= (l + z)(l,0,V p )> fc M = iu M , 1 + « = (3) 
where V p = dx p /dw are the transverse velocity components on the observer's sky. 




FIG. 2: Observational coordinates based on the past lightcones of the observer's worldline. 
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The metric of an arbitrary spacetime in observational coordinates is 

ds 2 = -A 2 dw 2 + 2Bdwdy + 2C P da; p dw + D 2 (dQ 2 + L PQ dx p dx Q ) (4) 

A 2 = (l + z)- 2 + 2C P V p + g PQ V p VQ, B = ^, (5) 

dy 

where v is a null afhne parameter, D is the angular diameter distance, and L P q determines the lensing distortion of 
images via the shear of lightrays, 

D 2 dL PQ 

The number of galaxies in a solid angle d£l and a null distance increment dv = Bdy is 

diV = Sn(l + z)D 2 Bdndy , (7) 

where S is the selection function and n — p m /m is the number density. 

Before specializing to isotropic observations, we identify how the observations in general and in principle determine 
the geometry of the past light cone w = wq of here-and-now, where y = z: 

• Given the intrinsic properties and evolution of sources, observations in principle determine: (a) the angular di- 
ameter distance D(wq, z, x p ), (b) the lensing distortion of images, L P q(wq, z, x r ). Thus the metric components 
g P Q on w = wq are determined. 

• Given (ll|) and the selection function, the number counts N(wo, z, x p ) in principle determine 
B(w , z,x p )p m (w ,z,x p ). 

• In principle, observations over extended timescales determine the instantaneous transverse velocities 
V p (wo, z,x Q) of discrete sources. Transverse velocities of clusters are in principle determined by the polar- 
ization of scattered CMB photons [8] . 

It follows that, in principle and for idealized observations: 

Idealized data {u^ , Bp m , #pq} on w = wq. (8) 

But this is insufficient to determine the geometry of the past lightcone, because we need go P = C p and we cannot 
separate out B and p m . Without gravitational field equations, we are unable to fully determine the spacetime and 
matter on the past lightcone, even assuming perfect information from discrete-source observations. As a consequence, 
it is also impossible to test gravity theories directly: 

Cosmological testing of gravity theories 

Even with perfect observations, we cannot determine the spacetime geometry and matter on our past 
lightcone without gravitational field equations. Thus observations cannot directly test GR on cosmolog- 
ical scales, or test any alternative theories of large-scale gravity, without making assumptions about the 
spacetime geometry. 

If GR holds, then ^ is exactly what is needed 0E]: 

Matter observations metric and matter on lightcone 

Observational data ^ is exactly the information needed for Einstein's equations to determine B and C p 
on w = wq, so that the metric and matter are fully determined on the lightcone. Then the past Cauchy 
development of this data determines g^, it M , p m in the interior of the past lightcone. 

If the matter observations are isotropic, then we can prove isotropy of spacetime and matter [SJIB]: 

Isotropy of lightcone matter distribution for one observer —> isotropic spacetime 

In an expanding dust region with A, if one fundamental observer measures isotropic angular diameter 
distances, number counts, bulk velocities, and lensing, 

then the spacetime and matter are isotropic, i.e. the region is Lemaitrc-Tolman-Bondi (LTB). 

It is an open question how this result translates to the case of almost-isotropy of observations, i.e. does this lead to 
almost-isotropy of the spacetime? 
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B. Isotropic CMB for one observer 

If the CMB is isotropic for one observer comoving with the matter, then along the observer's worldline C we have 



[see (171] 



f(x c ,p) = F{x c ,E), F M1 ... W L =0 = F fll ... fle L for l> 1. (10) 



In other words, all covariant multipoles of the distribution function beyond the monopole [see (A32)], and their time 



derivatives, must vanish along C. By (A33), the radiation momentum density (from the dipole) and anisotropic stress 
(from the quadrupole) vanish: q^\c = = tt^Ic- However, without the Copernican assumption, we are not able to 
deduce directly the vanishing of spatial derivatives V u F lll ... lle | c , and then we cannot show isotropy of the spacetime 
geometry about C. 

Isotropy of the CMB alone is not sufficient, since the matter could in principle be anisotropic, even if this is physically 
unnatural. In order to rule out artificial counter-examples, it would be sufficient to characterize the minimal matter 
isotropy that combines with CMB isotropy to give geometric isotropy. 

If we adopt the Copernican Principle, then the CMB alone leads to a powerful result, as discussed in the next 
section. 

III. WHAT IS THE BASIS FOR SPATIAL HOMOGENEITY? 

Homogeneity cannot be directly observed - we are effectively unable to move away in cosmic time or distance 
from here-and-now and hence cannot probe spatial variations on constant time slices; effectively, our observations 
only access the past lightcone of here-and-now. Direct observation thus cannot distinguish between an evolving 
homogeneous distribution of matter and inhomogeneity with a different time evolution - since the past light cone only 
accesses a 2-sphere in each constant-time slice (see Fig. [IJ. Thus we are forced to adopt the Copernican assumption. 
We first consider matter observations, then an exactly isotropic CMB, and finally the case of an almost-isotropic 
CMB. 

A. Isotropic matter observations for all observers 

If all observers see isotropy of the 4 matter observables, then we have geometric isotropy along all worldlines and 
thus homogeneity follows. 

Isotropy of lightcone matter distribution for all observers — > FLRW 

In an expanding dust region with A, if all fundamental observers measure isotropic angular diameter 
distances, number counts, bulk velocities, and lensing, then the spacetime is FLRW. 

This is an improved form of the Cosmological Principle - based on isotropy of specific observables and not on 
assumed geometric isotropy. 

In fact there is a much stronger statement than this, based only on one observable, the angular diameter distance, 
and only for small redshifts [3]: 

Isotropic distances up to 0(z 3 ) for all observers —> FLRW 

If all fundamental observers in an expanding spacetime region measure isotropic angular diameter distances 
up to third-order in a redshift series expansion, then the spacetime is FLRW in that region. 

A covariant proof of this IIP] is based on a series expansion in a general spacetime, using the method of Kristian 
and Sachs [11]. The redshift may be expanded in terms of the angular diameter distance: 

z = [PrV^^+^^^^V.V^J^ 2 
1 



G 



k ll k v k 'k f) V IJl VuV a Uf ) + -kfkrWkPRapVpUv] ^D 3 + 0(L> 4 ). (11) 



o 



Here all terms are evaluated at the observer O in the unit direction e M of the lightray, where [see pfc]: 

k» = -(i + z )( u M + e M) ) u^e" = 0, e^e" = 1. (12) 



() 



The covariant derivative of the 4- velocity is decomposed as in (A5|. The O(D) coefficient is the observed Hubble rate 
(generalizing the FLRW quantity), 



e 



Isotropy at lowest order, for all observers O, thus enforces vanishing acceleration and shear 



= = a^ v . 



Using (14 1, the next coefficient becomes 



1 



[k^k v k a V i y v u a } = - (26 2 + Pm - 2A) - 2V / ,6e" + (E^ + w^w^e'V 



(13) 



(14) 



(15) 



where V M is the covariant spatial derivati ve (|A3[ ) , E^ v is the electric part of the Weyl tensor ( A8 1 , and the angled 
brackets denote the spatial tracefree part (A2|. Isotropy at 0(D 2 ) therefore imposes homogeneity of the expansion 
rate, V M = 0, and the condition E UI , — —uji^uj u \. However, this condition is identically satisfied by virtue of (14) 



and the shear propagation equation (A14). The constraint equations (A18), (A19) then show that curlw^ = and 
Hpv — V^w^, where is the magnetic Weyl tensor. The 0(D 3 ) coefficient imposes further constraints, and we 



find that = V M /> 



E^v = H^. Putting everything together, we have a covariant characterization of FLRW. 



B. Exactly isotropic CMB for all observers 



It is commonly assumed that isotropy of the CMB for all observers leads obviously to FLRW, without the need for 
any proof. In fact, the proof is far from obvious, and requires a nonlinear analysis of the general Einstein-Liouville 
equations. The starting point is a pioneering mathematical result by Ehlers, Geren and Sachs [12] . They assumed that 
the only source of the gravitational field was the radiation, i.e. they neglected matter and A. This can be generalized 
to include self-gravitating matter and dark energy [I] (extending [13TfT8] ): 



CMB isotropy for all observers 

In a region, if 



FLRW 



• collisionless radiation is exactly isotropic, 

• the radiation four-velocity is geodesic and expanding ; 

• there is dust matter and dark energy in the form of A, quintessence or a perfect fluid, 
then the metric is FLRW in that region. 

The fundamental 4-velocity is the radiation 4-velocity, which has zero 4-acceleration and positive expansion: 



A„ 



o, e > o. 



(16) 



Isotropy of the radiation distribution about means that photon peculiar velocities are isotropic for comoving 
observers; thus in momentum space, the photon distribution depends only on components of the 4-momentum p 11 
along u M , i.e., on the photon energy E = —u^p^: 



f(x,p)=F(x,E), F w ... w = 0far*>l. 



(17) 



In other words, all covariant multipoles of the distribution function beyond the monopole, defined in (A32|, must 
vanish. In particular, (A33) shows that: 



(18) 



Equation (17) also implies that the radiation brightness octupole Tluva an d hexadecapole Tluvap are zero. These are 



source terms in the anisotropic stress evolution equation, which is the 1 = 1 case of (A37). The general nonlinear 
form of the ^ v evolution equation is O [20] 



3 15 5 
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u a ^ va& = 0. 



(19) 



7 



Isotropy removes all terms on the left except the third, and thus enforces a shear-free expansion of the fundamental 
congruence: 



0>«, = . 



(20) 



We can also show that u 11 is irrotational as follows. Together with (16 1, momentum conservation for radiation, i.e., 
(A30) with I = r, reduces to 



V M p r = . 



(21) 



Thus the radiation density is homogeneous relative to fundamental observers. Now we invoke the exact nonlinear 
identity for the covariant curl of the gradient, (A6|: 



curl Vupt = -2p r w„ => QprUin = . 



(22) 



where we have used the energy conservation equation (A29) for radiation. By assumption > 0, and hence we deduce 
that the vorticity must vanish: 



uj^ = 0. 



Then we see from the curl shear constraint equation (A19) that the magnetic Weyl tensor must vanish: 

H, lv = . 



(23) 



(24) 



Furthermore, (21) actually tells us that the expansion must also be homogeneous. From the radiation energy 
conservation equation (A29l, and using (18 1, we have — —Zp-[/Ap T . On taking a covariant spatial gradient and 
using the commutation relation (A7), we find 



v M e = o. 



(25) 



Then the shear divergence constraint, ( A18 ), enforces the vanishing of the total momentum density in the fundamental 
frame, 



<z" = o =► £t?(p;+pSK = o. 



(26) 



The second equality follows from (A26), using the fact that the baryons, CDM and dark energy (in the form of 
quintessence or a perfect fluid) have vanishing momentum density and anisotropic stress in their own frames, i.e., 

q*^ = = tt*^ , (27) 

where the asterisk denotes the intrinsic quantity (see Appendix A). If we include other species, such as neutrinos, 
then the same assumption applies to them. Except in artificial situations, it follows from ( 26 ) that 



< = 0, 



(28) 



i.e., the bulk peculiar velocities of matter and dark energy [and any other self-gravitating species satisfying (27)] are 
forced to vanish - all species must be comoving with the radiation. 

The comoving condition ( 28 ) then imposes the vanishing of the total anisotropic stress in the fundamental frame: 

^ : = 5>r = J2^(p*i+p>>i } = o. (29) 



where we used (A27), (27) and (28). Then the shear evolution equation (A14) leads to a vanishing electric Weyl 
tensor 



E„ u = 0. 



(30) 



Equatio ns (|2 6|) and (|29|), now lead via the total momentum conservation equation (A12) and the iS-divergence con- 

(31) 



straint (A20), to homogeneous total density and pressure: 

V> = = V^p . 
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Equations (16), (20 1, (24), (125]), (26), (129]) and (31) constitute a covariant characterization of an FLRW spacetime. 



This establishes the EGS result, generalized from the original to include self-gravitating matter and dark energy, and 
presented in a fully covariant form. It is straightforward to include other species such as neutrinos. The critical 
assumption needed for all species is the vanishing of the intrinsic momentum density and anisotropic stress, i.e., 
(27 1. Equivalently, the energy- momentum tensor for the /-component should have perfect fluid form in the /-frame. 



The isotropy of the radiation and the geodesic nature of its 4-velocity then enforce the vanishing of (bulk) peculiar 
velocities Vj . We emphasize that one does not need to assume that the matter or other species are comoving with 
the radiation - it follows from the assumptions on the radiation. 

In fact this result can be dramatically strengthened on the basis of a theorem by Ellis, Treciokas and Matravers 
[19]: we do not need vanishing of all multipoles, but only the dipole, quadrupole and octupole! The key step is to 
show that the shear vanishes, without having zero hcxadccapolc the quadrupole evolution equation (19) no longer 



automatically gives a^ v — 0, and we need to find another way to show this. The ETM trick is to return to the 



Liouvillc multipole equation (A36). The 1 = 2 multipolc of this equation, with Fn = F^ v = F^ va = 0, gives 



12 d 
63dE 



We integrate over E from to oo, and use the convergence property E 5 F flua ^ — > as E — > oo. This gives 

f°° dF 
<W„ E5 dE dE = °- 



(32) 



(33) 



Integrating by parts, the integral reduces to —5 J °° E 4 FdE. Since F > 0, the integral is strictly negative, and thus we 
arrive at vanishing shear, a^u = 0. Then the proof above proceeds as before, incorporating matter and dark energy 
to extend the ETM result [J: 

CMB partial isotropy for all observers — > FLRW 

In a region, if 

• collisionless radiation has vanishing dipole, quadrupole and octupole, 



F = F — F =0 



(34) 



• the radiation four-velocity is geodesic and expanding, 

• there is dust matter and dark energy in the form of A, quintessence or a perfect fluid, 
then the metric is FLRW in that region. 

This is the most powerful observational basis that we have for background homogeneity and thus an FLRW back- 
ground model. 



C. The real universe: almost-isotropic CMB 

In practice we can only observe approximate isotropy (see Fig. [3]). Is the previous result stable - i.e., does almost- 
isotropy of the CMB lead to an almost-FLRW Universe? This would be the realistic basis for a spatially homogeneous 
Universe (assuming the Copernican Principle). It was shown to be the case, but subject to further assumptions, by 
[13 [21]: 



A more realistic basis for homogeneity 

In a region of an expanding Universe with dust and cosmological constant, if all observers comoving with 
the matter measure an almost isotropic distribution of collisionless radiation, and if some of the time and 
spatial derivatives of the covariant multipoles are also small, then the region is almost FLRW. 

We emphasize that the perturbative assumptions are purely on the photon distribution, not the matter or the 
metric - and one has to prove that the matter and metric are then perturbatively close to FLRW. Once again, a 
nonperturbative analysis is essential, since we are trying to prove an almost-FLRW spacetime, and we cannot assume 
it a priori. 
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FIG. 3: The CMB temperature from WMAP 7-year data, showing anisotropies at the level ~ 10" 5 . (Credit: NASA/ WMAP 
Science Team.) 



Almost-isotropy of the photon distribution means that 

F Mv .. Mt (x,E)=0(e), 1>1, (35) 

where e is a (dimensionless) smallness parameter. The brightness multipoles TlM e have dimensions of energy density 
and we therefore normalize them to the monopolc II = p r /47r, so that HmJ^I = 0(e). 

The task is to show that the relevant kinematical, dynamical and curvature quantities, suitably non-dimensionalized, 
are 0(e). For example, the dimensionful kinematical quantities may be normalized by the expansion, er^/G, 
The proof then follows the same pattern as the proof above of the exact result - except that at each stage, we need 
to show that quantities are 0(e) rather than equal zero. 

However, in order to show this, we need smallness not just of the multipoles, but also of some of their derivatives. 
Smallness of the multipoles does not directly imply smallness of their derivatives, and we have to assume this [T5II22) . 
It remains a difficult open problem whether these additional assumptions may be removed. If all observers measure 
small multipoles, then it may be possible to use almost-isotropy of the matter distribution to show that the time and 
space derivatives on cosmologically significant scales must also be small. 

A number of experiments has been proposed to test the Copernican Principle by looking for violations of isotropy 
at events down our past lightcone, as discussed in the next section. The almost-isotropic CMB result then gives a 
framework for probing inhomogeneities via such observations. Indeed, these tests may provide a way of constraining 
spatial gradients of the low-£ multipoles. 

In addition, it may be possible to strengthen the almost-isotropic result above by proving that it is sufficient for 
only the first 3 multipoles and their derivatives to be small. This would represent a more realistic foundation for 
almost- homogeneity. 



IV. TESTING HOMOGENEITY 



Although we cannot directly probe homogeneity by observations, we can test for violations of homogeneity. If 
we find no violation, then the indirect evidence for homogeneity is strengthened. However, if even one significant 
violation is discovered, then homogeneity will have been disproved. There are two broad classes of observational tests 
of homogeneity, one based on consistency relations that hold in FLRW and the other based on using galaxy clusters 
as probes of the anisotropy seen in the CMB from positions at cosmological distances down our past lightcone. 



A. Consistency of distances and expansion rate 

The effective standard candles provided by supernovae observations lay the basis for a consistency test of homogene- 
ity. There are two geometric effects on distance measurements: the curvature bends null geodesies and the expansion 
changes radial distances. These are coupled in FLRW models via 
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and then one can combine the Hubble rate and distance data to find the curvature today: 

O _ HHz)[(l + z)D' L (z)-D L (z)?-(l + z) A 



This relation is independent of all other cosmological parameters, including dark energy - and it is also independent 
of the theory of gravity. It can be used at a single redshift to determine SIko- Furthermore, it is the basis for a CP 
test, proposed by |23j - since flxo is independent of z, we can differentiate ( |37[ ) to get the consistency relation: 

K{z) := (1 + zf + H 2 {z)[(l + zf {D L (z)D'l(z) D*(z)} + D 2 L (z)] 
+ (1 + z)H{z)H'{z)D L {z) [(1 + z)D' L {z) - D L {z)] 
= for FLRW geometry. (38) 

Note that JC — for any FLRW geometry, independent of curvature, dark energy, matter content, and theory of 
gravity. In realistic models we should expect |/C(z)| ~ 10~ 5 , reflecting perturbations from large-scale structure 
formation. Significantly larger values indicate a breakdown of homogeneity: 

JC(z) significantly different from => non-FLRW universe. (39) 

Carrying out this test should not be more difficult than carrying out dark energy measurements of w(z) from SNIa 
data, which require H'(z) from distance measurements or the second derivative D'l(z). 

This is the simplest test of homogeneity, and its implementation should be regarded as a high priority. Another 
test involves the time drift of the cosmological redshift 24J, but this will only be feasible on a much longer timescale. 



z+Az 




FIG. 4: Radial and transverse BAO scales. 



Finally the BAO feature itself provides in principle a homogeneity test. Future large-volume surveys will allow the 
detection of the BAO scale in both radial and transverse directions. The physical lengths in radial and transverse 
directions of a feature with redshift extent Az and subtending an angle A9 are 

where D is angular diameter distance and Hii is the expansion rate in the radial direction (see Fig. |4|. The two scales 
are equal for the BAO feature in an FLRW background, where expansion is isotropic at all points. Any significant 
disagreement between the radial and transverse BAO scales would signal a breakdown of remote isotropy and thus of 
homogeneity: 

i,, 

— - — 1 significantly different from => non-FLRW universe. (41) 
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B. Sunyaev-Zeldovich effect: temperature of scattered CMB photons 

Galaxy clusters with their hot ionized intra-cluster gas, act via scattering of CMB photons like giant mirrors that 
carry information about the last scattering surface seen by the cluster (see Fig. [5]). In other words, clusters allow us 
in principle to indirectly probe inside our past lightcone. 

CMB photons are scattered into our line of sight, thus inducing spectral distortions in the CMB temperature that 
we observe (25J [26]. The thermal Sunyaev-Zeldovich (SZ) effect from the thermal motion of electrons reflects the 
monopole seen by the cluster. If the blackbody temperatures of CMB photons arriving at the cluster from points 
inside our lightcone are significantly different from the blackbody temperature that we directly observe, then there 
will be a significant distortion. Such a signal would indicate that the cluster sees a significantly anisotropic CMB, 
hence violating the CP and homogeneity [27]. (This test has been applied to a class of LTB models by [28 .) 

The bulk radial motion of the cluster gas induces a kinetic SZ signal that reflects the CMB dipole seen by the 
cluster. If this is large, then there would be a violation of the CP and homogeneity. This has been applied to classes 
of LTB models by [251 EH] , but it is in fact more generally applicable as a test of homogeneity. 

In summary, 

Large (non-perturbative) thermal or kinetic SZ temperature effect non-FLRW universe. (42) 




FIG. 5: CMB photons from inside our past lightcone are scattered into our line of sight by the ionized gas in galaxy clusters, 
thus carrying information about the level of isotropy seen by the clusters. 



C. Sunyaev-Zeldovich effect: polarization of scattered CMB photons 

Analogous to the SZ effect on CMB temperature, there is an SZ effect on CMB polarization [5j. The cluster bulk 
transverse velocity, and the CMB monopole, quadrupole and octupole, as seen by the cluster, induce modifications 
in CMB polarization via scattering off the cluster gas. These effects have been computed in FLRW models [3TM33] . 
The SZ polarization signals potentially contain more information than the SZ temperature signals, and this has been 
proposed for the standard FLRW model as a way to map the CMB quadrupole seen at remote locations, thus lessening 
the cosmic variance [33H36] . 

From a more general viewpoint, the SZ effect on polarization is in principle a powerful probe of the CP and 
homogeneity In principle, if CMB polarization data indicate large transverse cluster velocities, or large modifications 
to the primordial polarization signal, this could signal a violation of the CP and homogeneity: 



Large (non-perturbative) SZ polarization effects =>■ non-FLRW universe. 



(43) 
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V. CONCLUSIONS 



It is important to re-examine the basic assumptions of the standard concordance model, especially in view of 
the problems raised by the fine-tuned and unnatural nature of dark energy. Here we have re-examined the central 
assumption of homogeneity. We argued that homogeneity cannot be established or directly confirmed by observations, 
given the inherent limitations of lightcone-based data. By contrast, isotropy can be directly probed by observations. 
We constructed the observational coordinates necessary to answer the question of which observables on the past 
lightcone are needed to prove isotropy. It turns out that 4 independent matter observables - angular diameter 
distances, number counts, lensing distortion and transverse velocities - are exactly sufficient to determine the spacetimc 
geometry. Isotropy of these 4 observables along one worldline imposes isotropy of the spacetime geometry. Surprisingly, 
isotropy of the CMB along one worldline does not in itself lead to isotropic geometry. 

To establish homogeneity, we are forced to adopt the Copernican Principle (CP). Isotropy of the 4 matter observables 
for all observers leads to homogeneity - giving an observational version of the Cosmological Principle. A surprising 
and more powerful result is that isotropy of only the angular diameter distances for all observers, and only to 0(z 3 ), 
enforces homogeneity. 

The strongest observational basis for homogeneity comes from the CP combined with the high isotropy of the 
CMB. The main additional assumption needed is that the CMB rest frame is geodesic. We outlined a covariant 
nonlinear proof and generalization of the original ECS result. And we highlighted the remarkable ETM result, i.e. 
that the vanishing of the dipolc, quadrupole and octupole is sufficient to enforce homogeneity. The realistic case, 
with almost-isotropy of the CMB, does not lead to almost- homogeneity without additional assumptions on some of 
the derivatives of the multipoles. It remains an open problem whether these assumptions can be avoided, possibly 
using further information from almost-isotropy of matter observations. 

Although we cannot directly observe homogeneity, we can test homogeneity, using observations that carry informa- 
tion from inside our past lightcone. We described how this can be achieved via consistency relations between distances 
and expansion rates, using supernova and BAO data, and via Sunyacv-Zeldovich effects from galaxy clusters on the 
CMB temperature and polarization. Up to now, none of these tests has yet indicated a breakdown of the CP and 
thus a violation of homogeneity But the further advance of high-precision data will provide new opportunities to 
apply and extend these critical tests. 

In summary: the standard homogeneous model of cosmology is successful, predictive and up to now robust against 
all the observational data, and against current tests of homogeneity The problem of a satisfactory explanation for 
dark energy may be resolved by advances in particle physics and quantum gravity. It is certainly reasonable to assume 
that the LCDM model is a good description of the Universe. Nevertheless it is also necessary to continuing probing 
the foundations of the model - not only the assumption of homogeneity, but also other critical questions, such as the 
problem of averaging and the need for a detailed understanding of light propagation in a lumpy Universe. 
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Appendix A: Nonlinear field and Boltzmann equations 

For convenience we repeat the summary of the 1+3 covariant Lagrangian formulation of the field equations and 
Boltzmann equation in full nonlinear generality, as given in [1]. For a chosen 4- velocity field the fundamental 
tensors are 

h[ii> Q^iv ~l~ ^^l^ui £fivot — Tliivafi^ \ (A-l) 

where projects into the instantaneous rest space of comoving observers, and e^ va is the projection of the spacetime 
alternating tensor ij^ap = — ^/—gS° 'i^S 1 v 8 2 a S 3 m. The projected symmetric tracefree (PSTF) parts of vectors and 
rank-2 tensors are 

V (ll) = hf/Vu , S ( ^) = {h (a a h„/ - \h^h^S a p . (A2) 

The skew part of a projected rank-2 tensor is spatially dual to the projected vector, — ^Sfj, ua S^ ua \ and then any 
projected rank-2 tensor has the decomposition — ^Sh^ + e Mt/Q 5 Q + S/^, where S = S a ph a ^ . The covariant 
derivative V M defines 1+3 covariant time and spatial derivatives: 

> = </"V„./" V,,./' . = //„ '/,"., • • • //,"V ../ ..... (A3) 

The projected derivative V M defines a covariant PSTF divergence, V M V^ , V^S^,,, and a covariant PSTF curl, 

curl = e^ c V v V a , curl = e a0{ ^ a S v) p . (A4) 

The relative motion of comoving observers is encoded in the PSTF kinematical quantities: the volume expansion 
rate, 4-acceleration, vorticity and shear, given respectively by 

9 = V M w M , An = iip, uj^ = curlw^, = V^u v ) V„?i M = ^Qh^ + e^ ua uj a + - A^u v . (A5) 

Key nonlinear identities are 



curlV^V := e^ <x V u V a ii = -2^uj ll , (A6) 
1 

-( 

3 



V( V -^) - V M ^> = V^m - ( o V + ^ + ^uk^ r ) V> ■ (A7) 



The PSTF dynamical quantities describe the sources of the gravitational field: the (total) energy density p — 
T^vU^u" and isotropic pressure p — \h llv T >lv (including A), momentum density = —T^ u u u , and anisotropic stress 
TT/j,!/ = T//xv) ; where is the total energy-momentum tensor. The Weyl tensor splits into PSTF gravito-electric and 
gravito-magnetic fields 

E» v = C m ^u a vP , H^ = ^e^C a ^ K u K . (A8) 

The Ricci and Bianchi identities, 

V[ M V„]U a = R^apu 13 , V^C^aff = — V^j-R^Q, — g-Rsvjaj, (A9) 

produce the fundamental evolution and constraint equations governing the covariant quantities. Einstein's equations 
are incorporated via the algebraic replacement R^ v = T^ v — \T a a g^ v . The resulting equations, in fully nonlinear 
form and for a general source of the gravitational field, are: 
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Evolution: 



p+(p + p)e + V% = -2i4% - a^ir^ , 

e + l& 2 + l(p + 3p) - VM M = -^a"" + + , 



09m + (P + P)A* + v mP + V" 71 ^ = + ^ V o^ v <t ~ A"np V , 

1 



-Swp + -curlA^ = ct^cj 



,/3 



1 



-i- 0-ff^,/ + cvAE^ - ^curl7i>„ = 3a a ^H v ) a - w Q £ Q(3(Al i/„)' 3 
3 1 



Constraint: 



2 ~ 



- curlw M - - V M + q^ = -2e, Ja/a u u A a 

curler^ + V^u„) - = -2A^u) v ) , 

!^ 1 

2 V ^-^+3 



1 . „ _„« 1 „ „« __ „ 1 



(A10) 
(All) 

(A12) 

(A13) 

(A14) 

(A15) 

(A16) 

(A17) 
(A18) 
(A19) 
(A20) 

(A21) 



The energy and momentum conservation equations are (A10) and (A12). The total dynamical quantities have 



contributions from all dynamically significant particle species: 

I" 4 " = T i V = P^u v + phf™ + 2q^u^ + if , Tf = pj« + pVt-T + 2q* (ll u I ) + iff , (A22) 
i 

where I = r,n,b,c, A labels the species. The asterisk on the dynamical quantities pj, ■ ■ ■ denotes that these quantities 
are measured, not in the u^-frame, but in the /-frame, whose 4-velocity is given by 

< = 7/ K + <) , = , 7/ = (1 - vjY 1 ' 2 . (A23) 

The fully nonlinear equations for the I dynamical quantities as measured in the fundamental u^-frame are: 

Pi = p* I + {ijv 2 I (p}+p}) + 2 1I q*/v I ^+7r* I ^v I ^v IlJ }, (A24) 

Pi = Pi + g [if v\ (p* +p}) + 2^iq*/vi fl + ir^vi^Viu j , (A25) 
Qi = qT + {p*i +P*i) v i + [ill - l )<lT ~ imTvivuL 1 + ~jiqTvi v v1 

+ 7r"? (P*i+P*i) < + **rv Iv ~ T$"*v Iv v Ict v»} , (A26) 
Trf = n*^ + { - 2u^7r*' J)a vi a + nr^vi^u" + 2^vfq* v) 

- 2jiqrvi a uK ,J I ) - ^* a0 vi a vi h^ + 7 ? (p*i + p*i) } . (A27) 

The terms in braces are the nonlinear corrections that vanish in the standard perturbed FLRW case. The total 
dynamical quantities in ( A10 1 |A21[ ) , are given by 



(A28) 
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Assuming that the species are non-interacting, they each separately obey the energy and momentum conservation 
equations lAlO} and (IA12I): 



Pi + (Pi + Pi)® + ^nQi = -ZA^qf - o^itf , 



1i 



iiu 



where the /-quantities are given by (A24|- |A27[ ). 

The covariant kinetic theory description starts by splitting the photon 4-momentum as 

p^ = E^ + e M ) , e% = 1 , e'X = . 



(A29) 
(A30) 



(A31) 



Here E = —u^p^ is the energy and e M = p^ / E is the direction, as measured by a comoving fundamental observer. 
Then the photon distribution function is decomposed into covariant harmonics via the expansion 



f(x,p) = f(x,E,e) = F + F^+F^e" + --- = J2 F M e (x,E)t 



(Mi) 



(A32) 



£>0 



where Me := [iifa ■ • • p>i and e Me := e Ml • • • e w . The PSTF multipoles Fm ( are a covariant alternative to the usual 
expansion in spherical harmonics. The energy-momentum tensor is T^ v {x) = J p^p v f(x,p)d 3 p, so the first 3 multipoles 
define the radiation dynamical quantities (in the u M -frame): 



p t = 4tt / E 3 F dE , q? = — / E 3 F f2 dE , <" = / E 3 F^ dE . 
Jo 3 J Q 15 J 

Thus n = (l/47r)p r ,IP = (3/47r)gf,II^ = (15/87r)7if where the brightness multipoles are 

IT,,,...,,, = / E F,., —a, dE . 



df 



The collisionlcss Boltzmann (or Liouville) equation is 



df 



df 



dv ^ dx^ a ^ ^ dp^ 



(A33) 



(A34) 



(A35) 



where p^ = dx^ /dv, and we neglect polarization. The covariant multipoles of df /dv are given by 



E\dv 



C + l) * 



(2£ + 3) 



(2* + 3) 

/>, ,f_ P a (l+l)(l + 2) pi— (£+2) rpl+3p I'^-AM 7 

^ e va{lil F Ml _ ±) - {2e + 3){2e + 5) E [E F^ Mi \ a 



2f. 



(2£ + 3)' 



E 1/2 E 3/2 F„ {AIe i <7^ e) u -E l 1 [E 2 e F {Me _ 2 ] ct w _ iW >, 



(A36) 



where a prime denotes d/dE. This is a fully nonlinear expression. Multiplying by E 3 and integrating over all energies 
leads to the brightness multipole evolution equations: 



o = n 



-en 



M, 



V (w n M,-i) 



(2^ + 3) 



(2£ + 3) 

„ (l-l)(l+l)(l + 2) 
(2£ + 3)(2J? + 5) 



V UuM, 



+ ' " ' ~ 2 U"iw, + (£ + ^A^n^^) 



(7" Q n 



volMi 



{21 + 3) 



(A37) 



The monopole equation is just the energy conservation equation, i.e., (A29) with I — r, the dipole is the momentum 
conservation equation (A30l, with I = r, and the quadrupole is (19). 



16 



[1] C. Clarkson and R. Maartens, Class. Quant. Grav. 2 7 (2010) 124008 [arXiv:1 005.2165 . 
[2] G. F. R. Ellis, Class. Quant. Grav. 28 (2011) 164001 [arXiv:1103.2335| . 

[3] P. Sarkar, J. Yadav, B. Pandey and S. Bharadwaj, Mo n Not. Roy. Astro n. Soc. 399 (2009) L128 |arXiv:0906.3431 
[4] F. Sylos Labini, Class. Quant.' Grav. 28 (2011) 164003 [arXiv: 1103. 5974] , 
[5] R. Maartens, PhD thesis, University of Cape Town (1980). 

[6] G. F. R. Ellis, S. D. Nel, R. Maartens, W. R. Stoeger and A. P. Whitman, Phys. Rept. 124 (1985) 315. 

[7] R. Maartens and D. R. Matravers, Class. Quantum Grav. 11 (1994) 2693. 

[8] R. A. Sunyaev and Y. B. Zeldovich, Mon. Not. Roy. Astron. Soc. 190 (1980) 413. 

[9] W. Hasse and V. Perlick, Class. Quantum Grav. 16 (1999) 2559. 

[10] C. A. Clarkson, Ph.D. thesis, University of Glasgow (1999) |[arXiv:astro-ph /00 08089 . 

[11] J. Kristian and R. K. Sachs, Astrophys. J. 143 (1966) 379. 

[12] J. Ehlers, P. Geren and R. K. Sachs, J. Math. Phys. 9 (1968) 1344. 

[13] W. R. Stoeger, R. Maartens and G. F. R. Ellis, Astrophys. J. 443 (1995) 1. 

[14] R. Treciokas and G. F. R. Ellis, Commun. Math. Phys. 23 (1971) 1. 

[15] J. J. Ferrando, J. A. Morales and M. Portilla, Phys. Rev. D 46 (1999) 578; 

[16] C. A. Clarkson and R. K. Barrett, Class. Quantum Grav. 16 (1999) 3 781 [arXiv:g r-qc/9906097 . 
[17] C. A. Clarkson and A. A. Coley, Class. Qua nt. Grav. 18 (200 1) 1305 |arXiv:astro-ph/0102172'P 
[18] S. Rasanen, Phys. Rev. D 79 (2009) 123522 [arXiv:0903.3013| . 

[19] G. F. R. Ellis, R. Treciokas and D. R Matravers, Ann. Phys. 150 (1983) 487. 

[20] R. Maartens, T. Gebbie and G. F. R. Ellis, Phys. Rev. D 59 (1999) 083506 arXiv:astro-ph/ 9808163 . 
[21] R. Maartens, G. F. R. Ellis and W. R. Stoeger, Phys. Rev. D 51 (1995) 1525 |arXiv:astro-ph/9501016| . 
[22] S. Rasanen, JCAP 0902 (2009) 011 |arXiv:0812.2872] . 

[23] C. Clarkson, B. Bassett and T. H.-C. Lu, Phys. Rev. Lett. 101 (2008) 011301 



arXiv:0712.3457 



[24] J.-P. Uzan, C. Clarkson and G.F.R. Ellis, Phys. Rev. Lett. 100 (2008) 191303 
[25] R. A. Sunyaev and Y. B. Zeldovich, Astrophys. Space Sci. 7 (1970) 3. 
[26] R. A. Sunyaev and Y. B. Zeldovich, Comments Astrophys. Space Phys. 4 (1972) 173. 

[27] J. Goodman, Phys. Rev. D 52 (1995) 1821 |arXiv:astro-ph/9506068] . 

[28] R. R. Caldwell and A. Stebbins, Phys. Rev. Lett. 100 (2008) 191302 [arXiv:071 1.3459] 
[29] J. Garcia-Bellido and T. Haugb0lle, JCAP 0804 (2008) 003 |arXiv:0807.1326|. 
[30] P. Zhang and A. Stebbins, Phys. Rev. Lett. 107 (2011) 041301 arXiv: 1009.3967 



arXiv:0801.0068 



[31] S. Y. Sazonov and R. A. Sunyaev, Mon. Not. Roy. Astron. Soc. 310 (1999) 765 arXiv:astro-ph/9903287 . 

[32] A. Challinor, M. Ford and A. Lasenby, Mon. Not. Roy. Astron. Soc. 312 (2000) 159 arXiv:astro-ph/9905227 . 

[33] M. Kamionkowski and A. Loeb, Phys. Rev. D 56 (1997) 4511 [arXiv:astro -ph/9703118|. 

[34] N. Seto and M. Sasaki, Phys. Rev. D 62 (2000) 123004 [arXiv:astro-ph/0009222|. 

[35] D. Baumann and A. Cooray, New Astron. R ev. 47 (2003) 839 |arXiv:astro-ph/0304416| . 

[36] E. F. Bunn, Phys. Rev. D 73 (2006) 123517 |arXiv:astro-p h/060327'TT7^ 



